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ABSTRACT 

We probe the possible anisotropy of the Universe by using the JLA compilation of 
type-la supernovae. We apply the Markov Chain Monte Carlo (MCMC) method to 
constrain the amplitude and direction of anisotropy in three cosmological models. For 
the dipole-modulated ACDM model, the anisotropic amplitude is consistent with zero 
at 68% C.L., and has an upper bound Ap> < 1.98 x 10 -3 at 95% C.L. Regardless of 
much larger uncertainty, we find the dipole direction of JLA is amazingly opposite to 
that of Union2. Similar results are found for the dipole-modulated u>CDM and CPL 
models. Thus, the Universe is still well consistent with the isotropy according to the 
JLA compilation. 
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1 INTRODUCTION 

The cosmological principle, which is one of the foun¬ 
dations in modern cosmology, says that the Univ erse is ho¬ 
moge neous and isotropic at large enough scales JPodelsonl 
1200311 . It is well consistent with the present observational 
data, such as the cosmic microwave background (CMB) 
radiation from Wilkinson Microwave Anisotropy P robe 


(WMAP1 (Bennett et al. 

20l| 

Hinshaw et al. 2013i') and 

Planck satellite (Ade et al 

2015a 

bj). Until now, the cosmo- 


logical observations are still in accordance with the cosmo¬ 
logical constant plus cold dark matter (ACDM) model which 
is based on the cosmological principle. Thus, the ACDM 
model becomes the leading model in modern cosmology. 

Despite the great successes it achieved, th e ACDM 
mode l still faces certain challenges llPerivolaropoulod 
l200Sl . l201ll ; iMariano fe Perivolaropoulosl 120131 1. As the 
improvements of accuracy, it is found from a large amount 
of observations that the Universe might deviate from 
statistical isotropy. These include the alignment of low 
multipoles in the angular power spectrum o f CMB tern per- 
ature fluctuations (iTegmark. de Oliveira-Costa fe Hamilton 


)ower asymm e try o f CMB temperat ur e anisotropy 
iBennett et al.l 1 20131 : lAde et akl 120141 : lEriksen et al.l 


2003 

: Bielewicz 

201C 

; Frommert 


12004 ICopi et al 
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20041: Hansen. Bandav fe Gorskil 120041 : lAkrami et al.l 


20 14 iQuartin fe NotarH 2015l l. the spatial variation of 


the el ectromagne t ic fin e -structure constant dWebb et al.l 
l201ll : iKing et al.1 120121 : iMolaro et al.l l2013h . the large- 
scale alignment of the quasar polarization ve ctors 
llHutsemekers fe Lamvl I200II : iHutsemekers et al.ll2005h . the 
large-scale bulk flow beyond the prediction of ACDM model 
d Kashlinskv et aPl2009l . 120101 : 1 Watkins. Feldman fe Hudsonl 
l2009lh and so on. All of these phenomena arouse us to 
rethink the validity of the cosmological principle. If the 
cosmological principle is proven to be failed, the modern 
cosmology should be rewritten. 

Due to their consistent absolute magnitudes, type-la 
supernovae (SNe la) are regarded as the ideal distance 
indicators to trace the accelerated expansion of the Uni¬ 
verse. In fact, they have been widely used to search for the 
anisotropic signals in the Universe. Especially, a statistic 
based on the extreme value theory s hows that the gol d data 
set is consistent wit h the isotropy ( Gupta. Saim_&Laskatj 
l2008h . The study dBlomqvist. Enander fe Mortselll 2010h 
on the angular covariance function of supernova magnitude 
fluctuations is consistent with zero dar k energy flu ctua¬ 
t ions by using the Union2 compilation (lAmanullah et al.l 
l2010ll . A “residual” statistic shows that the isotropic 
ACDM model is not consi stent with the U nion2 data 
with z < 0.05 at 2 — 3<r dColin et al.l l20lIL There are 
no significant evidence for deviations f rom th e isotropy 
in th e anisotropic Bianchi-I cosmo lo gy IlCampanelli et al.l 
l201ll : ISchucker. Tilouin fe Valenti gOldf) , Bi anchi-III 

and Kantowski-Sachs metrics llKoivisto et alJ 1201 il l. 
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and Randers-Finsler cosmolog y dChang et al.l l2014al . 
120131 : IChang. Li fe Wand l2013h . The hemisphere com¬ 
parison is used to study the Union2 data and shows 
certain preferred directions (jSchwarz_&^einhors^ 


Antoniou&Periv^aro£Oulo^__]2010|; [Cai&Tuc 


2007; 


2012 ; 


Kalus et al] 120131 : lYang. Wang fe Chul |20 l~dT ~ Chang fe Linl 
20151) . By dividing the Union2 supernovae into 12 subsets 


according to their galactic coordinates, a dipole of the 
deceleration parame t er is preferred at more than 2u level 
llZhao. Wu fe Zhangl 120131 ). By combining the data of 
Union2 and gamma -ray burst s, the isotropic ACDM model 
is well permitted (ICai et al.l 120131 ) while th e anisotropic 
Finsle r cosmology is preferred at arou nd 2cr dCha ng et alJ 
l2014bl ). By using the data of Union2.1 d Suzuki et al. I 20121 ) 
and gamma-ray bursts, a model-ind ependent way shows a 
dipolar anisotropy at more than 2er (IWang fe Wangll2014l) . 
It has been found that there may be certain correlation 
betwee n the fine structure dipole and the dar k ene rgy 
dipole dMariano fe Perivolaropoulosll20f^ : lLi et ~aI1l2015l ). A 
fully-Bayesian method was developed to remove the system- 
atics in the Union datasets a nd the anisotropic cosmology 
does not seem to be reflected dHeneka. Marra fc Amendolal 
120141) . 


Recently, a new sample of SNe la was released 
by the SDSS collaboration, which is call ed the “joint 
light- curve analysis” (JLA) compilation dBetouIe et al.1 
2014) . Compared to pr evious compilat ions such as Union 2 
l Amanullah et alJl2010h and Union2.1 dSuzuki et alJl20l3) . 
the number of SNe la in the JLA compilation is highly en¬ 
larged and the systematic uncertainties are significantly re¬ 
duced. Recently, the JLA SNe la have been used to probe 
the anisotropic Hu bble diagr am in Bianchi type I cosmology 
dSchucker. Tilqu in fe Valenti 1201 41) and test the cosmolog¬ 
ical princ iple dBengalv. Bernui fe Alcaniz I 2 OI 5 I ). However, 
the work dBengalv. Bernui fe Alcanidl2015 ) did not consider 
the full covariance matrix between SNe la. I 11 this paper, we 
use the JLA compilation to restudy the anisotropic Hub¬ 
ble diagram of the Universe. Unlike certain previous works 
which have neglected the correlations between any two SNe 
la, we make use of the full covariance matrix to construct the 
likelihood (or chi-square). In addition, we use the method of 
MCMC sampling in our analysis. It has been shown that the 
statistical significance of the previously claimed evidence for 
a preferred direction could be highly lowered if the full co- 
variance matrix of SNe la is considered in the Union2 com¬ 
pilation dJimenez. Salzano fe Lazkozll2015l ). We want to see 
whether the anisotropic signals in the accelerated expansion 
of the Universe still exist in the newly released JLA compi¬ 
lation. 

The rest of the paper is arranged as follows. In section 
0 we briefly introduce the JLA dataset, and present the 
anisotropic cosmological models and the numerical method 
used in our analysis. In section [3] we give constraints on 
the anisotropic amplitudes and directions for the anisotropic 
expansion of the Universe. Finally, our conclusions are given 
in section [4] 


2 DATA AND METHODOLOGY 

The anisotropic expansion of the Universe can be in¬ 
duced by assuming that the dark energy has anisotropic 


repul s ive force dArmendariz-Pi conl 120041 ; iKoivisto fe Motal 
|200SI: ISalehi fe Aftabil l2015i) . or the backgrou nd space- 
time has a certain pre f erred direction dChanget_aL| 
20131: IChang, Li fe Wane! 120131: iLi et aid 120131. l2015t 
Schucker, Tilquin fe Valentll2014l ). and so on. In this paper, 


we assume a dipole modulation to describe the Universe de¬ 
viating from the isotropic background. Phenomenologically, 
the direction-dependent distance modulus can be given as 


AHh = Mth (1 + A D (h ■ p)). 


(1) 


where Ad denotes the amplitude of the dipole modula¬ 
tion, h is the dipole direction, p is the unit 3-vector point¬ 
ing towards the supernova, and p,th denotes the theoret¬ 
ical distance modulus predicted by the isotropic ACDM, 
wCDM or CPL models. Here the anisotropic amplitude 
Ad is assumed to be a constant over the whole redshift 
range. In the galactic coordinates, the dipole direction n 
can be parameterized as (l,b), where l and b are the lon¬ 
gitude and latitude, respectively. In such a parametriza- 
tion, we have n = cos(b) cos(Z)i + cos(6) sin(Z)j + sin(6)k, 
where i, j, k are the unit vectors along the axes of a 
Cartesian coordinates system. The position of the ith su¬ 
pernova with galactic coordinates (li,bi ) can be written as 
Pi = cos (bi) cos(Zi)i + cos(fei) sin(Z;)j + sin(6;)k. 

In the spatially-flat isotropic background spacetime, we 
can express the luminosity distance dz, (z) of a supernova in 
terms of the redshift z, 


d L {z) = 


1 + z 
Ho 


F 

Jo 


dz' 

W) 


( 2 ) 


where Hq = lOO/i km s -1 Mpc -1 is the Hubble constant, 
and E(z) is a function of redshift. The isotropic distance 
modulus /fth can be given by 


/Xth = 51 Ogl 0 lfe' (3) 

In equation (f2|l. the quantity E(z) is a function of redshift 
a. The expression of E(z) depends on a specific cosmological 
model. I 11 the ACDM model, it can be expressed as 


E 2 (z) — Hm(l + z ) 3 + (1 — Sim) , (4) 


where is the energy density of matter today. In the 
wCDM model, it can be expressed as 

e 2 (z) = n m (i + z) 3 + (i - n m )(i + z ) 3(1+w) , ( 5 ) 

where w = p/p denotes the equation of state of dark en¬ 
ergy. I n the Chevallier-Polarski- Li nder (C PL) parametriza- 
tion llChevallier fe Polarskill200ll : lLindedl2003l) , the equation 
of state of dark energy is redshift-dependent, and it is pa¬ 
rameterized by w = wo +wiz/(1 + z). In this case, E(z) can 
be expressed as 

E 2 (z) = n m (l+ 2 ) 3 + (l-n m )(l+ 2 ) 3(1+, " 0+, " l) exp f—3wi 

\ 1 + Z 

(6) 

In this p aper, we use the most recently published JLA 
compilation dBetoule et all [2014) of SNe la to constrain 
the anisotropy of the Universe. The JLA compilation con¬ 
sists of 740 well-calibrated SNe la in the redshift range of 
2 € [0.01, 1.30]. It is a collection of several low-redshift sam¬ 
ples, all three seasons from the SDSS-II, three years from 
SNLS, and a few high-redshift samples from the Hubble 
Space Telescope (HST). All of the SNe la have high-quality 
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light curves, so their distance moduli can be abstracted with 
high precision. The positions of SNe la in the sky of equato¬ 
rial coordinates system can be found at the website of IAU 
Central Bureau for Astronomical Telegram^. To compare 
with others’ work, we transform the positions of SNe la into 
the galactic coordinates. 

From the observational point of view, the distance mod¬ 
ulus of a SN la can be abstracted from its light cur ve through 
the empirical linear relation (IBetouIe et al.M 2014f 

A = m *B ~ {Mb — a x Xi + /3 x C), (7) 


where m* B is the observed peak magnitude in rest frame B 
band, Mb is the absolute magnitude depending on the host 
galaxy properties complexly, Xi is the time stretching of 
the light curve, and C is the supernova color at maximum 
brightness. The three light-curve parameters m B , X\ and 
C are different from one supernova to other one and can 
be derived directly from the light curves. The two nuisance 
parameters a and /3 are assumed to be constants for all the 
supernovae. 

For the JLA samples, the isotropic luminosity distance 
of a supernova can be given as 


d /. (^hel ■ Zcmb) 


1 + «hel f Zcmb dz' 

Ho Jo E(z') ’ 


( 8 ) 


where z c mb and Zhei denote the CMB frame redshift and 
heliocentric redshift, respectively. Then we can obtain the 
anisotropic distance modulus fit h in equation ©• Using the 
observed distance modulus A in equation 0 , the anisotropic 
cosmological models can be fitted to the JLA dataset by 
using the chi-square as 

Xjla = (A - Mth) 1 C” 1 (A - Mth), (9) 


where C is the covar iance matrix of A, and it is presented in 
IBetouIe et al.l (120141 ). 

In order to directly compare with pr evious works, we 
also apply our method to the Union2 (lAmanullah et all 
l20inh dataset. The Union2 data set consists of 557 SNe la 
with well-observed redshift in the range of z € [0.015, 1.4], 
The distance moduli and their uncertainties are extracted 
from the SALT2 light-curve fitter. The directions of SNe la 
are w ell loc alized in th e sky of the equatorial coordinates 
llBlomavist. Enander fe Mortsellll2010l) . The data in LTnion2 
are usually assumed to be uncorrelated. Thus, the chi-square 
of Union2 is simplified to 


XUnion2 f ) 


fl — flth 


( 10 ) 


where the summation runs over all the SNe la data. 

In this paper, we employ the Markov Chain Monte 
Carlo (MCMC) method to estimate the model parameters. 
The joint likelihood is given by £ ex exp(—X 2 /2). The nui¬ 
sance parameters such as a and /3 are marginalized. We mod¬ 
ify the publicly available Cosmologic al Monte Carlo sampler 
(CosmoMC) (iLewis fe Bridle! 12002) to estimate the back¬ 
ground parameters and anisotropic parameters. For ACDM 
model, the isotropic parameter can be well constrained by 
supernovae data. For ruCDM and CPL models, however, the 
isotropic parameters can’t be well constrained by using the 


1 http://www.cbat.eps.harvard.edu/lists/Supernovae.html 


parameters 

ACDM 

i«CDM 

CPL 

Ad 

i{°] 

6[°] 

< 1.98 x 10 -3 

Q1 fi+ 107 
c>iD -110 
c+41 
°-60 

< 2.09 x 10~ 3 

320~*” 107 

ozu_ 104 

_4+45 

4 -61 

< 2.05 x 10~ 3 

3is±m 

0+36 

°-54 


Table 1. The 95% upper bound of dipole amplitude Ad, and the 
preferred direction (l, b) with la uncertainty in three cosmological 
models. 


supern ovae data only. Following the method in lBetoule et al.l 
(t2014l) . therefore, we combine the supernovae data with 
the Planck 2013 results of CMB temperature anisotropy 
(Planck2013), the WMAP9 observations of CMB polariza¬ 
tions (WP), and the SDSS-III BOSS DR11 Baryon acoustic 
oscillations data (BAO) to constrain the isotropic param¬ 
eters. Once the isotropic parameters are given, we can fix 
them and fit the anisotropic parameters with JLA SNe la 
only. 


3 RESULTS 


As was mentioned above, we study the anisotropic sig¬ 
nals of the dipole-modulated ACDM, wjCDM and CPL mod¬ 
els by using the JLA sample. The nuisance parameters such 
as a and f3 can be marginalized, since they are not model 
parameters with significant meanings. We just focus on 
studying the anisotropic signals, thus neglect the topic of 
model comparison. Our final results are listed in Table [T] 
where the models and the anisotropic amplitudes and pre¬ 
ferred directions are given. The likelihood distributions of 
the anisotropic parameters Ad, l and b in three cosmological 
models are plotted in Figured] In the last panel of Figured] 
we also plot the distribution of Xjla- 

In the modulated ACDM model, the isotropic param¬ 
eter can be well constrained by using the JLA d ataset 
only. The result is 12 m = 0.295 ± 0.034 (IBetouIe et al.l 
[2014! ). which is computed for a fixed fiducial value of Ho = 
70 km s" 1 Mpc~ 3 . The change of Ho does not affect the 
best-fit value of 12 m• By using the MCMC approach, and 
fixing 12 m at 0.295, the anisotropic amplitude is constrained 
as Ad < 1.98 x 10~ 3 at 95% C.L., which is consistent with 
the isotropy within la uncertainty. This implies that the 
JLA compilation shows no significant evidence for the devi¬ 
ations from isotropy. In the galactic coordinates, the dipole 
direction points towards F (l, b ) = (316° , —5° igpo) at 

68% C.L. By contrast, the Union2 sample gives constraints 
on the dark energy dipole as Ad = (1.3 ± 0.6) x 10“ 3 and 
{l, b) = (309.4° ± 18.0°, -15.1° ± 11.5°) at 68% C.L., which 
implies that the anisotropi c expansio n of the U niverse is per - 
mitted at more than 2<r (lMariano fc Perivolaropoulosll2O0 ) . 
However, the statistical significance can be highly lowered 
if the full covariance matrix o f SNe la is consi d ered in the 
Union2 compilation d Jimenez. Salzano fe Lazkodl2015h . Our 
constraints on the anisotropic parameters are listed in the 
second column in Table [T] Their likelihood distributions are 
illustrated by the solid curves in Figure [T] 

In the modulated wCDM model, there is an extra 
background parameter, i.e., the equation of state of dark 
energy w. This parameter can’t be well constrained by 
using the JLA data only. The combined constraint from 
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Figure 1 . Likelihood distributions for the amplitude Ad and 
direction (/, b) of dipole modulation in three cosmological models. 
The distribution of Xjla a ^ so showed in the last panel. 


Planck2013+WP+BAO+JLA gives Q, M = 0.303 ± 0.012, 
■ w = -1.027 ± 0.055 and H 0 = 68.50 ± 1.27 km s -1 Mpc -1 
dBetoule et alJ 1201 40. By fixing these parameters at their 
best-fitting central values, we use the MCMC approach to 
constrain the anisotropic parameters. Similarly to the above 
discussions, we obtain the 95% upper bound on anisotropic 
amplitude Ad < 2.09 x 10 —3 , and preferred direction (l, b) = 
(320° 1jo4° ’ -6i°) 68% C.L., which are listed in the 

third column in Table [l] Their likelihood distributions are 
illustrated by the dashed curves in Figure [T] 

In the modulated CPL parametrization, there are 
two extra background parameters, i.e., the equation of 
state of dark energy parameters wo and wi. These 
two parameters also can’t be well constrained by us¬ 
ing the JLA data only. The combined constraint from 
Planck2013+WP+BAO+JLA gives = 0.304 ± 0.012, 
w 0 = -0.957 ± 0.124, wi = -0.336 ± 0.552 an d H 0 = 
68.59 ± 1.27 km s -1 Mpc -1 dBetoule et al.ll2014T ). By fix¬ 
ing these parameters at their best-fitting central values, the 
anisotropic amplitude is constrained to be Ad < 2.05 x 
10 -3 at 95% C.L., and the preferred direction (l, b) = 
(318° , —8° 1 3 4 o) at 68% C.L. These results are listed 

in the last column in Table Q] Their likelihood distributions 
are illustrated by the dash-dotted curves in Figure [lj 

From Table |T| we can see that all the three cos¬ 
mological models give consistent results. At first glimpse, 
the dipole dire ctions we o bt ained seem to be consis¬ 
tent with that of lMariano fc Perivolaropoulosl d2012lf . How¬ 
ever, this is in fact not the case. In our calculation, 
we constrain dipole amplitude to be non-negative, i.e., 


Ad ^ 0, and let the dipole dire c tion r uns over the whole 
sky. iMariano fc Perivolaropoulosl (|2012h parameterized the 
dipole modulation as fit h = Ath(l — Ao(n ■ p)), which has 
a sign difference from our parametrization of equation ©. 
Therefore, o ur dipole directio n is actuall y oppo site to that 
obtained bv IMariano fc Perivolaropoulosl (j2012l h 

To test if the discrepancy between our results and previ¬ 
ous works is due to the different dataset or different method. 
We apply our method to Union2 dataset, such that di¬ 
rect comparison with previous results can be made. The 
constraint on isotropic ACDM model from Union2 date- 
set gives Qm = 0.274 ± 0.040, with fiducial parameter 
Ho = 70 km s -1 Mpc -1 . Then we use the MCMC ap¬ 
proach to constrain the anisotropic parameters. We obtain 
A d = (0.54t°; 13 ) x 10 -3 and (l,b) = ( 142 °+"°, 11 °+^°) 
at 68% C.L. Regardless of the much larger uncertainty, 
th e dipole direction we obtai n ed is consistent with that 
of lMariano fe Perivolaropoulosl ll2012l h However, our results 
show that the Union2 dataset is consistent with isotropy 
at la C.L., and the 95% upper bound on anisotropy is 
Ad < 1-40 x 10 -3 . On the other hand, if we apply the 
least-square method, we obtain Ad = 1.10 x 10 -3 , and 
( l,b ) = (126°, 18°). Both the anisotropic amplitude and 
preferred direction are well consis tent with the results of 
IMariano fc Perivolaropoulosl ll2012l h 


4 CONCLUSIONS 

In this paper, we probed the possibly anisotropic ex¬ 
pansion of the Universe by using the recently released 
JLA compilation of SNe la. We considered the dipole- 
modulated deviation from the isotropy in three different 
dark energy models. We obtained similar constraints on the 
anisotropic amplitude and direction in three cases. This in¬ 
dicates that the preferred direction anisotropy is insensitive 
to the isotropic background models. Especially, our MCMC 
studies show that the anisotropic amplitude has an upper 
bound D < 1.98 x 10 -3 at 95% C.L., and the dipole di¬ 
rection points towards (l, b) = (316° 4gJo) for the 

dipole-modulated ACDM model. These results imply that 
the there is no significant evidence for anisotropy in the JLA 
dataset. For comparison, we also applied MCMC method 
to the Union2 dataset, and we got Ad < 1-40 x 10 -3 at 
95% C.L. The dipole direction of the Union2 points towards 
( l , b) = (142° , 11° ^ 27 °) a t 68 % C.L., which is consistent 

with previous results. We surprisingly found that the dipole 
direction derived from the JLA is approximately opposite to 
that from the Union2. 
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